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With the use of the Hubbard model bound hole states in YBa 2 Cu 3 6+2 ,. are 
studied. For the parameters of this crystal the exchange interaction between 
the spin-carrying chain ion 0~ and Cu — O plane sites is shown to ensure 
the formation of a large ferromagnetically ordered clusters around holes in 
the plane. 



1. Introduction 



Over the past decade the magnetic structure of cuprate perovskites has been 
studying very intensively in the hope that it might provide insight in the 
physical origin of high-temperature superconductivity. It is presently well 
established that Cu0 2 planes of undoped and lightly doped cuprates are an- 
tiferromagnetically ordered. Some experiments in such crystals, for example 
magnetization measurements in YBa 2 Cu30 6+:r [HJ, suggest the existence of 
local ferromagnetically ordered regions around holes on the antiferromagnetic 
background of Cu0 2 planes. In application to another type of crystals such 
excitations were first considered in Ref. 2 and were called ferrons. Ferron 
models of high-T c materials were discussed in Refs. 3,4. The explanation 
of the ferron formation, suggested in Ref. 1, is based on the comparison of 
antiferromagnetic and ferromagnetic interactions created by a hole localized 
on an oxygen orbital. This idea goes back to work [0] where it was used for 
the explanation of the observed destruction of the long-range antiferromag- 
netic order at small hole concentrations. The mechanism suggested in the 
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latter work is essentially based on the supposition that nonhybridized oxy- 
gen and copper orbitals are appropriate starting states for the consideration. 
However as shown in Ref. 6, low-lying states are strongly hybridized. As a 
consequence, the mechanism of Ref. 5 is not appropriate for explaining the 
destruction of the long-range ordering. For similar reasons this mechanism 
may not be responsible for the ferron formation. 

Notice that the mentioned strong hybridization of oxygen and copper 
orbitals and the crystal structure of CuC>2 planes open the way to reduce 
the Emery Hamiltonian for these planes to the more simple t-J or Hubbard 
Hamiltonians with only one effective orbital per plane sell ||. In the t-J 
model the ferron formation was considered both by the exact diagonalization 
of small clusters [| and in an infinite plane with the use of the spin-wave 
approximation [|H], [UJ . Both approaches showed that in the periodic lattice 
ferrons are energetically unfavorable for parameters of cuprate perovskites. 
The next step is to check whether point defect can stabilize ferrons. In our 
recent paper fl2| we have considered local states generated by the interstitial 
2 ~ ion in L^CuCU+s and found that this type of defects is unable to stabi- 
lize ferrons. In YE^CusOe+z, where ferrons were experimentally observed, 
for small x oxygen ions start to fill chain sites and can be considered as in- 
terstitial defects above the CuC>2 plane (see Fig. 1). Their valency is equal 
to -2 or -1 [O and in the latter case the defect has the spin 1/2. It can 



be expected that, in contrast to the spinless O ion, such type of defects 
can stabilize states with a large spin, i.e. ferrons, due to an arising exchange 
interaction. 



2. Basic Formulae 



To check this idea we use the one-band Hubbard model for the description 
of a CuC>2 plane. The model Hamiltonian reads 

H = UJ2 n V n U ~ t J2 c L c i+a<x, (1) 

1 laa 

where the operator c\ a creates an electron in the effective orbital on the 
plane site 1 (the orbitals are centered at copper sites) with the spin a, a = 
(0, ±a), (±a, 0), where a is the in-plane copper distance, n\ a = c\ a ci a , U and 
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t are the Coulomb repulsion and the hopping matrix element, respectively. 
A chain oxygen 0~ is situated at 01 position and we add the following terms 
to the Hamiltonian 

H 1 = a l a a ~J2tm (a£<W + c m a a o) + U p a\a^a\ai, (2) 

a ma 

where aj. is the p-electron creation operator in the 0~ chain site, A is the 
difference in chain and plane site energies, U p is the Coulomb repulsion on 
the chain oxygen, and t m is the hopping matrix element between chain and 
plane sites. Summation over m proceeds over the plane sites in the vicinity 
of the "defect" chain oxygen site (these plane sites are shown in Fig. 1) and 
we suppose that t m « t for these sites. 

For YBa2Cu306+z the parameters in Eqs. (1) and (2) can be estimated 



as t m « t « leV, A « -2eV, U « 10eV,U p « 4.5 eV TO, III. For such 



parameters the correlation terms in Eqs. (1),(2) and the first term in Eq. (2) 
can be used as the unperturbed Hamiltonian and the second, hopping terms, 
as the perturbation in the perturbation theory expansion. The resulting 
effective Hamiltonian reads 

H cS = P ( -t ^2 c\ a ci +aa + t;Yj s i+aSi + J m S m S J P, 

V lacr ^ la m / 

U> m " {A + U P )(U-AY [ ' 

where 

are the components of the spin operators Si on the plane sites. The analogous 



Hamiltonian was derived in [15]. The components of the spin operator S Q 
on the chain site are defined analogously through the operators a a . P is the 
projection operator on the subspace where the chain level is singly occupied 
and the plane levels are singly occupied or empty. The first two terms in 
Eq. (3) represent the one-band t-J Hamiltonian on a square lattice. The 
summation over 1 proceeds over all plane sites, while the summation over m, 
over the "defect" plane sites. In Eq. (3), the "defect" term can be represented 
as the sum of terms with S z and S ± spin components. The term with S z 
components plays the main role in the ferron stabilization, whereas the terms 
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with S ± components lead only to additional lowering of energy levels. For 
this reason for simplicity we take into account only components with S z in 
the "defect" term in Eq. (3). Let us suppose that the 01 site is occupied 
by an electron with S z = 1/2. In this case the effective Hamiltonian can be 
rewritten in the form 



H. 



J 



off 



lacr 



la 
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QZ 
Or, 



(4) 



Since the sign of J m is the same for all sites in the "defect" region, the last 
term can stabilize ferrons. 



We use the spin-wave approximation |T0|, [L6| for the description of spin 
excitations in the plane. The approximation reduces to neglecting terms of 
the third and higher orders in the spin-wave operators 6 m introduced into 
Eq. (4) by the Holstein-Primakoff formulas 



-b\b, 



1/2 



where $i = (l - frfo) ' , II = (n/a, it /a), and Pf 
Let us introduce the hole creation operator h\ 



[1 + cr exp(im)]/2. 

E CT P\C\ a for the Neel 
state \M) . This classical Neel state can be determined by the condition 
h\J\f) = (for the second Neel state analogous formulas are obtained by 
substituting Pf with 1 — P\). After the unitary transformation which di- 
agonalizes the Heisenberg part of the Hamiltonian [JTUJ (the second term in 
Eq. (4)) the effective Hamiltonian acquires the form 



ll'a 



H eS = £ t [jiiW4&i-i'(wi'+a + vv) + H.c. + ^ £ oj v b\b\ +v 

1 ii' 



- £ (iimim + J2 m b\ n b m ) exp(illm), 



(5) 



where u\, v\ and uj\ are the Fourier transforms of cosh(a k ), — sinh(a k ), and 

= 4yl - t£, respectively, a k = (1/8) In [(1 + j k )/ (1 - 7 k )] , 7 k = 
[cos (k x a) + cos (k y a))/2, n m = h^hm. Due to the fast decrease of uy +a + Vy 
and uj\i with the growth of |1'| only the components with 1' = (±a, 0), (0, ±a) 
for the sum in the first, kinetic energy term and the components with 1' = 
(0, 0), (±o, ±a), (±2a, 0), (0, ±2a) for ojy in the second, magnetic energy term 
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are retained in Eq. (5) in the subsequent calculations. In this equation, J lm = 
Jm/4, J 2m = 2(ul + vl)J lm where u 2 + vg « 1.16. 

To investigate the lowest eigenvalues and eigenfunctions of Hamiltonian 



(5) we use the modified Lanczos algorithm [[K], [T^J. In one step of this 
algorithm we determine a final state |f) from an initial state |i) according to 
the procedure 

(i|i) = i, Si = (i\n\i), v\i) = (n-s l )\i). (6) 

Here V is the normalization constant which is determined from the condition 
(f |f) = 1. After one Lanczos step the energy is minimized in the subspace 
of the states |i) and |f) (it follows from (6) that (i|f) = 0). The obtained 
function = ci|i) + C2\f) which minimizes the energy is used as the initial 
function in the next Lanczos step. This procedure is continued until the 
necessary accuracy of the eigenvalue is reached. If the initial function |i) is 
characterized by some spin projection S z , |f) also corresponds to the same 
S z . Therefore, it is possible to obtain the lowest eigenvalues of Hamiltonian 
(8) with a given S z by starting from a state with the respective S z . Notice 
that states with S z > 1/2 contain flipped spins near the hole and correspond 
to ferrons. 

The number of spin configurations which is necessary for describing the 
states |i) and |f) grows from step to step. To avoid the overflow of the com- 
puter storage and to reduce the computation time, in each Lanczos step the 
state is restricted to the j configurations with the largest amplitudes []10| 



After normalization this restricted state is used as the initial state for the 
next recursion step of the Lanczos procedure (6). We used j = 200. The 
|f) state obtained from such state and used for the calculation of the 
eigenvalue contains of the order of 10 4 spin configurations. This is compara- 
ble to the number of configurations used in exact diagonalization studies of 
small lattices. To check the accuracy of the obtained results in some cases 
we increased j up to 250 which led only to small quantitative changes (of the 
order of 10 _4 t) in the eigenvalue. 
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3. Results and Discussion 



Using the above- described algorithm we have calculated the lowest one-hole 
states characterized by S z < 13/2. For parameters of cuprate perovskites the 
ratio J/t can be estimated to lie in the range 0.2 - 0.5 fllQf . According to 



Eq. (3) the value of J m /t for nearest neighbours of 01 oxygen site can be 
estimated to be equal 1. For other sites in the "defect" region the values 
of J m /t = 0.7 is accepted. The results of our calculations for the lowest 
one-hole eigenvalues are shown in Fig. 2. As seen, the interaction between 
the 0~ ion spin with plane spins leads to the ferron formation: the lowest 
state has S z = 7/2. The optimal ferron size can be expected to grow for two 
and larger number of adjacent 0~ ions in a chain. 

In summary, with the use of the spin-wave approximation and the Lanczos 
algorithm the localized hole states induced by an 0~ chain ion in YBa2Cu306 +x 
were investigated. The lowest hole state is characterized by the z component 
of the total spin S z = 7/2 which corresponds to a comparatively large ferro- 
magnetic region with 3 overturned spins around a hole. This means that the 
spin-carrying 0~ ion can stabilize ferrons. 
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Figure captions. 

Fig. 1. The "defect" region in the vicinity of the chain 0~ ion (01). 
Fig. 2. The dependence of energies of the localized hole states induced 
by the 0~ ion on S z . 
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